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ABSTRACT 

We employ a sequence of Doppler images obtained with the Coronal Multi-channel Polarimeter 
(CoMP) instrument to perform time-distance seismology of the solar corona. We construct the first 
k — uo diagrams of the region. These allow us to separate outward and inward propagating waves and 
estimate the spatial variation of the plane-of-sky projected phase speed, and the relative amount of 
outward and inward directed wave power. The disparity between outward and inward wave power 
and the slope of the observed power law spectrum indicate that low-frequency Alfvenic motions suffer 
significant attenuation as they propagate, consistent with isotropic MHD turbulence. 
Subject headings: Sun: solar wind - Sun: magnetic fields - Sun: corona 
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1. INTRODUCTION 

The Sun's corona is the primary source of the dis- 
turbances to the Earth's environment known as space 
weather. Despite its relevance, our knowledge about fun- 
damental physical parameters in the solar corona is still 
quite rudimentary. The measurement of coronal tem- 
peratures, densities and abundances is difficult. And al- 
though magnetism dominates the force balance in the 
corona, there exist very few direct measurements of the 
strength and orientation of coronal magnetic fie lds (see, 
e.R.,tB rosius & White 200l iLin. Kuhn fc Cou lter 2004; 
iTomczvk et al. 2008). Consequently, the mechanisms re- 
sponsible for heating the corona, driving the solar wind, 
and initiating coronal mass ejections remain poorly un- 
derstood. 

Over the past decade a great deal of progress has 
been made in observing and interpreting the signa- 
tures of magneto-hydrodynamic (MHD) waves in the 
corona. Various types of MHD waves have been identi- 
fied including fas t mode shocks seen as EIT waves (e.g., 
[Thompson et al.l 119981 ), slow magneto-acoustic (MA) 
waves in polar plun ies (iDeForest fc Gurm an 1998) and 
coronal loops (e.g.. i De Moortel. Ireland fc Walsh 2000: 
Marsh fc Wal sh 2006) a nd fast MA kink wave s (e.g., 
Aschwanden et al.. ,19991 : iNakariakov et all fT999h . The 



last of these has been the focus of much recent work 
where Transition Region and Coronal Explorer (TRACE; 
[Handy et al. 1999) measurements of impulsively excited 
transverse coronal loop displacements have been used 
to constrain the strength of the coronal mag netic field 
and o ther atmospheric parameters (p \.schwande n et al.l 
11999'; 'Nakariakov et al.' '1999'; 'Aschwanden et al." "2002|; 
ISchrijver, Aschwanden, Title 2002). Such investigations 
of the coronal magnetic structure and plasma environ- 
ment, in the limit of MHD, has been dubbed coronal seis- 
molo g y (e.g.. I Asch wanden 2002i: [Nakariakov fc Verwichtd 
120051 : iBan erjee et al. 2007). Observable loop displace- 
ment oscillations are relatively rare and only a few dozen 
loops have been analyzed seismologically. 
Recently, we discovered spatially and temporally ubiq- 
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uitous waves in the solar corona in the time series of 
Doppler images taken with the Coron al Multi-channe l 
Polarimeter (CoMP) instrument (Tomc zvk et al.l [20071 ). 
These data present the exciting possibility to greatly 
widen the application of coronal seismology. Analysis re- 
vealed propagating waves with an amplitude of ~300 m/s 
rms (1 km/s peak-to-peak) and peak power at a period 
of 5 minutes with Alfvenic phase speeds (~1 Mm/s) and 
with very small associated intensity fluctuations. These 
waves were interpreted bv Tomczvk et al.. (2007.) as prop- 
agating Alfven waves. 

These ubiquitous propagating waves eluded earlier in- 
vestigators for several reasons. First, these waves do not 
have appreciable intensity fluctuations so they have not 
been seen with intensity imaging instruments. Also, a 
velocity amplitude of 1 km/s at a period of 5 minutes 
corresponds to a loop displacement amplitude of 48 km 
which is nearly an order of magnitude smaller than the 
365 km (0.5'') pixel size of TRACE and nearly two or- 
ders of magnitude sm aller than the ^2000 km pix el size 
of SOHO/EIT (2.6''; iDelaboudiniere et al.l[T995[ ). The 
observation of these oscillations as loop displacements is 
probably beyond the capability of current instruments. 
In addition, Doppler imaging measurements with the re- 
quired sensitivity level, spatial extent and temporal ca- 
dence are difficult (or impossible) to obtain with ex- 
isting slit scanning spectrograph instrun ients. For ex- 
ampl e, the SOHO/SUMER instrument (|Wilhelm et all 
Il995f ) has a spectral sampling of approximately 10 km/s 
(43 mA) per pixel in first order; observation of these 
ubiquitous coronal waves with SUMER would require 
line centroiding to considerably better than 0.1 pixel. 
Recent analysi s of s pectra from the EIS instrument 
(|Culhane et al.l [2007) on Hinode CKosugi et al] 120071 ) 
show a portion of the corona oscillating at a period 
very nea r 5 minutes with a pea k-to-peak amplitude of 
2 km/s (|Van Doorsselaere et al.l [200 8). indicating that 
these ubiquitous waves may be marginally observable 
with EIS. Ground-based spectroscopic observations with 
a coronagraph showed clear evidence of oscillations, but 
their identificatio n as propagating Al fvenic disturbances 
was inconclusive (|Sakurai et al.|[2QQ2[ ). 
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Fig. 1. The Fe XIII 10747 A intensity (panel A) and derived 
wave propagation angles (panel B) in the region of the corona stud- 
ied by CoMP on 2005 October 30. On panel A we show, for later 
reference, the active loop region (dashed trapezoid) over which we 
have summed Fourier power spectra of the Doppler velocity to 
construct Fig. [2] the sample path used to determine the wave diag- 
nostics (dashed curved line;) and the surrounding (dotted) region 
to examine the spatially averaged k — uj diagram studied in Fig. 0] 
The wave propagation angles are measured counter-clockwise from 
due East. 

In this paper we exten d the analysis of th e same 
dataset used in the work bv lTomczvk et alJ (|2QQ7( ). They 
noted the clear presence of outward and inward propagat- 
ing Alfven waves in the corona, but no attempt was made 
to separate them. We present here a method to separate 
outward and inward propagating waves through the con- 
struction and analysis of coronal time series projected 
onto the paths of wave propagation. These time-distance 
diagrams are suitable for spatio-temporal Fourier analy- 
sis allowing us to build the first k — uo diagnostic diagrams 
of the corona. As we will show, these allow significant 
improvement in the determination of the phase speed of 
the waves, which enhances our ability to perform coronal 
seismology. 

In the following section we will briefly discuss the 
CoMP instrument and the observations analyzed in this 
paper. In Sect. 3 we discuss the data analysis techniques 
that we employ to derive the time-distance diagnostic 
measures. In Sect. 4 we discuss the implications of our 
results and place them in context with previous work. 
We conclude with a discussion of the potential for future 
observations to advance the field of coronal seismology. 

2. INSTRUMENT AND OBSERVATIONS 

The CoMP instrument is a combination polarimeter 
and narrowband tunable filter that can measure the com- 
plete polarization state in the vicinity of the Fe XIII 
infra-red coronal emission lines at 10747 A and 10798 
A and the chromospheric 10830 A He I line with a 2.8 
Rq field-of-view and a sampling of 4.5 arcseconds/pixel. 
The bandpass of the filter is tuned electro-optically with 
a bandwidth of 1.3 A (FWHM). CoMP was deployed be- 
hind the 20-cm aperture Coro nal One Shot coronagraph 
(|Smartt, Dunn fc FisherllTosTh on the equatorial spar of 
the Hilltop facility at the Sacramento Peak Observatory 
of the National Solar Observatory in Sunspot New Mex- 
ico. 



As in iTomczvk et al.l (|2007[ ), the observations used in 
this paper were obtained on Oct 30, 2005 and were re- 
stricted to the linear polarization (Stokes I, Q and U 
only) of the 10747 A Fe XIII line, sampling three wave- 
lengths (10745.2, 10746.5, and 10747.8 A) across the line. 
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Fig. 2. — The log-log Fourier spectrum of the measured Doppler 
velocity averaged over the region shown in panel A of Fig. ^ Also 
shown are the low- (red), high-frequency (blue) power-law fits and 
the photospheric p-mode spectrum (green). 

The exposure time for the images was 250 ms and image 
groups of full polarization and wavelength information 
were obtained at a cadence of approximately 28.7 sec- 
onds. The analysis presented in this paper is restricted 
to the first 350 image sets due to their superior quality, 
covering 2.8 hours starting at 14.261 UT. We selected a 
sub-array on the east limb for further analysis (see. Panel 
A of Fig. [T]) which contains both active region loops and 
a coronal cavity. 

For each pixel in the selected sub-array we fit a Gaus- 
sian profile to the intensity vs. wavelength and obtained 
estimates for the central intensity, central wavelength, 
and linewidth. The line-of-sight (LOS) component of the 
velocity was derived from the Doppler shift of the central 
wavelength, and we computed the plane-of-sky direction 
of the magnetic field (0) using the measured Stokes Q 
and U at line center and the relationship 
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In Fig. [2] we show a power spectrum of the coronal 
LOS velocity fluctuations averaged over the area above 
the active region to the South of the equator in Fig. [H It 
shows that the power spectrum has a peak at a frequency 
of 3.2 mHz; above 3.2 mHz the signal is dominated by os- 
cillatory power. The spectrum of photospheric p-modes 
computed by averaging all intermediate degree p-mode 
spectra with spherical harmonic degrees between and 
199 observed on one day with the LOWL instrument 
(|Tomczvk et al.l ri995f ) is also shown on the plot^. We 
notice that both the low- and high-frequency portions of 
the spectrum obey power-law scaling. Fitting linear rela- 
tionships to the log-log coronal wave power spectrum in 
the low- (0.1 > f > 2.5 mHz) and high-frequency (3 > f 
> 10 mHz) ranges we obtain values of the slopes of -1.54 
± 0.02 and -1.49 ± 0.02, respectively, where we quote 
the 1-cr errors in the fit gradients. The relevance of the 
power-law scaling in the wave spectrum will be discussed 
below (see. Sect. [5]). 

3. ANALYSIS 1. DETERMINATION OF WAVE 
PROPAGATION DIRECTION 

The first step of the time-distance analysis is to 
determine the direction of wave propagation. The 

^ The p-mode spectrum shown has been filtered to remove power 
below 1 mHz and been scaled to fit on the plot. 
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analysis method used here is describe d in detail in 
iMcIntosh. De Pontieu TomczvS (|2QQ8D and is a mod- 
ificati on of the method employed in iTomczvk et alJ 
(|2QQ7l ) . We determine the coherence (the frequency space 
equivalent of the cross-correlation) of the velocity time 
series of each pixel with its surrounding pixels. These 
coherence maps reveal an "island" of high coherence (de- 
fined as coherence greater than 0.5). We define the di- 
rection of wave propagation to be aligned with the is- 
land and compute the orientation of the island to be 
along a line which minimizes the sum of perpendicu- 
lar distance from the po ints in that island to the line. 
The method employed in iTomczvk et al.l (|2QQ7l ) differed 
in that it determined the time domain cross-correlation 
and fit the orientation of the islands using a simple linear 
least-squares fit. We subsequently determined that the 
original method failed to determine the wave angle well 
when the wave angle was highly inclined. 

The map of wave angles is shown in Panel B of 
Fig. [1] along with the coronal Fe XIII intensity. Com- 
parin^ this an^ l e map with the corresponding panel of 
ITomczvk et al.l (|2QQ7l ) (Fig. 4 panel B; showing angles 
with errors of less than 5°) we see that while the general 
spatial pattern is the same, the new approach has been 
able to reach considerably higher inclinations of the wave 
trajectory. 

In Figure [3] we show the comparison of the wave direc- 
tion with the directly measured magnetic field direction 
from the linear polarization measurements and Eq. [H 
Due to the above algorithm refinement, the correlation 
between the wave direction and magnet ic field orienta - 
tion i s improved over that presented in ITomczvk et al.l 
([2QQ7h . The magnetic field direction obtained from linear 
polarization measurement s is subject to the well know n 
90° Van Vleck ambiguity (IVan Vlecklll925l : lHouselll972h . 
This ambiguity occurs when the angle between the in- 
cident radiation field and the magnetic field direction 
passes through an angle of 54.6°. The ambiguity is ev- 
ident in the figure as the locus of points that lie par- 
allel to but 90° below the correlation line. We define 
the points which are subject to the Van Vleck ambigu- 
ity as those that are more than 45° away from the cor- 
relation line and find that about 20% of the points in 
the figure are subject to the Van Vleck ambiguity. Ne- 
glecting these points, the correlation coefficient between 
the wave propagation direction and the magnetic field 
direction is 0.91. Th is confirms our previous assertion 
(|Tomczyk et al.l 120071 ) that these waves are propagating 
along magnetic field lines imbedded in the corona. We 
anticipate that the Van Vleck effect will, in the future, 
provide a strong constraint on the magnetic field geom- 
etry and the location of the emitting region along the 
line-of-sight. 

4. ANALYSIS 2. CONSTRUCTION OF SPACE-TIME AND 
K -oj DIAGRAMS 

We employ the map of wave directions computed above 
to trace the wave paths in the corona, starting from any 
location in the field of view. This is accomplished by 
stepping a small amount (typically one pixel) from a 
starting point in the direction of wave propagation and 
then interpolating a new direction at that point and re- 
peating the process. An example wave path is shown in 
Fig. □ with a length of 80 pixels (--250 Mm). Next, we 
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Fig. 3. — Scatter plot of the measured magnetic field direction 
(panel F of Fig. 1 in Tomczvk et al. 2007) and the derived wave 
propagation angles (panel B of Fig. [T}. The solid red line indi- 
cates one-to-one correspondence while the blue line designates the 
material subject to the Van Vleck ambiguity that is offset by 90°. 



compute the velocity for each point along the wave path 
for each velocity image in the time series and stack them 
into a space-time diagram. We use a cubic interpolation 
algorithm to map the velocities in the original data cube 
onto the wave path. An example space-time diagram is 
shown in panel A of Fig. [4] for the wave-path structure 
traced in Fig. [H We see outward and inward propagat- 
ing waves as linear structures in the time series that are 
inclined to the right and left, respectively. The angle of 
the structure to the vertical in the space-time diagram 
gives a measure of the wave's phase speed. 

The k — uo diagram of the path is the Fourier power 
spectrum of the space-time diagram, see Panel B of 
Fig. m It shows ridges of power in both positive and 
negative uo that correspond to outward and inward prop- 
agation of the waves, respectively. We can also see that 
the wave power in each direction is not equal. In this case 
the outward directed wave power is about a factor of 2 
greater than the inward power. The ridges of power in 
positive and negative uo appear to have the same inclina- 
tion, i.e., the same phase speed. The dashed lines shown 
in the panel correspond to phase speeds of 600 km/s 
whi ch is significantly low er than the 1.2 Mm/s reported 
by (|Tomczvk et al.l 120071 ) for this region. Also visible in 
the panel is a vertical ridge of power at very low values of 
uo. This is due to the large amount of low frequency ve- 
locity power consistent with the /"^^^ spectrum shown 
in Fig 2. The example single-path k — uo diagram can 
be enhanced by averaging the diagrams of neighboring 
wave-paths. In Panel C of Fig. [H we have averaged the 
k — uo diagrams for the paths in the dashed region sur- 
rounding the single path shown in Fig. [TJ We see that 
this region-averaged k — uo diagram shows the same prop- 
erties as that of Panel B of Fig. [H 

Fitting the ridges in these k—uo diagrams to extract the 
phase speed of the waves in each direction is not straight- 
forward due to the large amount of low spatial and tem- 
poral frequency power present, especially for the single- 
path timeseries, and so we must find another method. 
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Fig. 4. — CoMP space-time diagram for the sample wave path 
traced out in Fig. [T] (panel A), the k — uo diagram for that space- 
time diagram (panel B) and the average k — uj diagram for the 
region around the sample wave path (panel C). 

Taking the simplest possible approach, we compute the 
space-time diagram of outward and inward waves by 
masking the positive and negative frequency halves of the 
complex k — ix) transform separately and compute the in- 
verse transform. The outward and inward filtered space- 
time diagrams corresponding to the example of Fig. [H 
are shown in panels B and C of Fig. O The outward and 
inward timeseries show a few key properties. First, the 
correspondence between the wave patterns in the filtered 
(panels B and C) and the unfiltered (original, panel A) 
timeseries is clear with the vertical stripes oriented as 
we would expect, with the outward waves inclined to the 
right and the inward waves inclined to the left. The 
phase speed of the waves, indicated by the gradient of 
the stripes, appears to be approximately constant over 
the duration of the time series in each case. 

Using the filtered timeseries we are able to measure 
the wave phase speed along the path by cross-correlating 
the timeseries at the center of the path with those across 
the entire length of the path. We determine the lag of 
the cross-correlation functions [where cc > 0.5] with a 
parabola (the peak of the parabola determines the lag 
of the wave) and obtain the phase speed by fitting the 
slope of the observed lag with position along the path. 
Examples of the cross-correlation functions and resulting 
space-lag plots for the outward, inward and original un- 
filtered timeseries are shown in Fig. [6l We see that the 
phase speeds of the waves in the outward, inward, and 
unfiltered timeseries are 587 (± 2), 596 (± 20) and 728 
(±6) km/s, respectively. We note that the outward and 
inward phase speeds are consistent, but the phase speed 
of the composite is much higher. This discrepancy is 
understandable in terms of linear wave theory by consid- 
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A) Sample Track Space-Time Diagram 
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B) Outward Wave Timeseries 
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Fig. 5. — Sample CoMP space-time diagram for the sample 
structure traced out in Fig. [T] (panel A), the corresponding k — 
UJ filtered outward (negative-o;; panel B) and inward (positive-o;; 
panel C) space-time diagrams. 
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Fig. 6. — Cross-correlation phase-speed determination for the 
reconstructed (outward - panel A; inward - panel B) and the un- 
filtered (panel C) space-time diagrams. The phase-speed in each 
case is the gradient of the linear-fit to the wave lag at different 
positions along the wave path. 



ering a situation where waves are traveling in opposite 
directions on the same structure with the same speed 
but with different amplitudes. When the amplitudes are 
equal the combination of the two traveling waves is a 
standing wave; as the relative amplitude increases in one 
direction or the other the net wave will appear to travel 
in the direction of dominant amplitude. Since the work 
of iTomczvk et all (|2QQ7[ ) did not separate outward and 
inward waves but instead analyzed the composite time 
series, that work severely overestimated (by a factor of 
2) the phase speeds except where the wave power was 
dominated by either the outward or inward propagat- 
ing com ponent - the ph ase speeds shown here supersede 
those in Tomczvk et al. (2007). 

As in Tomczyk et al. (2007) we can repeat all the steps 
of the analysis pixel by pixel over the field-of-view shown 
in Fig. [TJ In doing this, we have computed wave paths 
of 31 pixels in length (~100 Mm) centered on the pixel 
of interest. We chose to use a short path for this step 
of the analysis to keep the wave paths from extending 
appreciably beyond the field-of-view. The value of 31 
pixels was the shortest length that gave reliable results. 
From these wave paths, we compute filtered space- time 
diagrams and extracted the phase speed in each, as well 
as that in the composite timeseries. We also compute the 
power in the outward and inward portions of the k — co 
diagram for each path. Results of this analysis are pre- 
sented in Figs. and [HI The panels of Fig. [71 show the 
resulting maps of outward and inward wave power as well 
as their ratio [ p°^^~p^^ ] from left to right while the pan- 
els of Fig. [8] show the respective phase speed maps in the 
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three cases (outward - panel A; inward - panel B; unfil- 
tered space-time diagram - panel C) and their respective 
errors. The power maps show excess power at the highest 
parts of the field-of-view. This is due to the decreasing 
coronal intensity at these heights and correspondingly 
higher noise levels there. The power ratio map shows 
that the corona is dominated by outward propagating 
waves except in a small region on the east side of the coro- 
nal cavity where inward wave propagation slightly dom- 
inates. Due to the dominance of outward wave power, 
the inward phase speed is generally less well determined 
than the outward phase speed and that is indicated by 
the larger errors seen. 

5. DISCUSSION 

The observed peak of the spectrum of coronal ve- 
locity perturbations in Fig. [T] shows an exact coinci- 
dence with the peak of the p-mode power spectrum. 
From this we conclude that the ubiquitous coronal 
waves observed in this frequency range are driven by 
solar p- modes. Recent work supports the contention 
that p-modes can be efficiently transported into the 
corona via coupling with magneti c fields rooted in the 
photosphere (De Pontieu, Erdelv i fc De Moortell l2QQ5l : 
iDe Moortell^ Rosneni2QQ71 : lHindman fc Jainll2QQ8^ 

We observe a power-law spectrum in the Doppler 
velocity over our entire frequency range, with excess 
power in the high-frequency (f > 2 mHz) oscillatory re- 
gion of the spectrum. A power-law spectrum has been 
widely observed in fluctuations of the solar wind with 
a slope between -5/3 and -3/2 which has been inter- 
preted as evidence of a turbulent cascade (see, e.g., 
iGoldstein RobertsllToOOh . A slope of -3/2 i s predicted 
for is otropic incompressible MHD turbulence (|Kraichnanl 
Il965[ ). The CoMP observations reveal a slope of -3/2 in 
both the low and high frequency regimes consistent with 
a turbulent cascade. 

It is noteworthy that the -3/2 slope of the high fre- 
quency (f > 3 mHz) tail of the coronal velocity spectrum 
is in contrast to the photospheric velocity spectrum in 
this frequency range which declines as ~ f~^'^ • This is 
evidence that the high frequency coronal spectrum is not 
merely imprinted by the photospheric spectrum. Either 
some mechanism modifies the photospheric spectrum as 
waves propagate into the corona or the coronal spectrum 
is modified in situ. 

We observe wave phase speeds that are approximately 
constant over the duration of the observations. This 
implies that the magnetic structures supporting the 
waves and the driving mechanism are likewise stable 
over this time scale. We note the large differences 
between inward and outward wave power, even along 
closed loop structures like the coronal cavity region. 
This can only be possible if the waves are damped 
appreciably in situ. The approximate time scale for 
damping must be of the order of the travel time along 
the loop, which amounts to just a few wave periods. 
This is consistent with the observed damping times 
of coronal loop oscillations ( Aschwanden et aiT [2QQ2). 
When considering the efficiency of these waves in heating 
the corona, or driving the solar wind (De Pontieu et al. 
[200 7l ). we must explore this damping process further, 
consider the origins of the observed large non-thermal 
line widths in the corona and exploit the developing 



self-consistent theories which involve the effective 
dissipation of low-frequency Alfy e nic motions (e.g;., 
"Cranmer & van Ballesooiien '2005"; 'Suzuki fc Inutsukal 
2005; Cranmer, van Ballegooiien fc EdgaJ l20Q7l : 
Verdini fc Vellil 120071 ). This wih be the subject of 
a forthcoming paper. 

As we have noted above, iTomczvk et al.l (|2007l ) 
interpreted the waves seen in the CoMP data as 
Alfven waves. This conclusion is consistent with 
the observations that the measured phase speeds of 
the waves are much larger than the sound speed, 
that the waves propagate along field lines and that 
the wa ves exhibit ver y small intensity p e rturba tions. 
Van Doo rsselaere. Nak ariakov fc Verwichtd (|2008l ) and 
Hindman fc Ja in (2008) suggest that these waves are 
more likely fast MA kink modes. The classical Alfven 
wave formulation (| Alfven [ [1942) assumes a uniform mag- 
netic field in a uniform density environment. MA kink 
modes arise in a mathematical formulation that assumes 
a cylindrical loop geometry with a den sity jump between 
the inside and outside of the loop (jEdwin fc RobertsI 
Il983l ). In such a formulation, the Alfven mode is tor- 
sional and will not be detected in observations where the 
loops are not spatially resolved, which is the case for 
the CoMP observations. In fact, coronal loops have yet 
to be conclusively resolved by any instr ument. We ad- 
mit the possibility that the waves seen in lTomczyk et al] 
(|2007) and this paper may, in fact, be fast MA kink 
modes. However, we have no information about the den- 
sity structure of the coronal features studied here, so an 
unambiguous interpretation in terms of a specific geom- 
etry is unwarranted at this time. Characterization of 
these waves as the more generic Alfven waves is in our 
opinion somewhat less presumptuous. 

Further, in the MA kink mode interpretation, under 
the usual assumption of a uniform magnetic field inside 
and outside the loop, the phase speed, Vphase^ is approx- 
imately given by the kink speed, Vkinki and is related to 
the Alfven speed inside the loop, Vai (= 
relation 



yphase ^ Vkink — ^Aii ' — \ 1 — 7 (2) 

y 1 + po/pi 

where pi and po are the densities inside and outside the 
loop, respectively and ii j is the mean molecula r weight 
inside the loop (see, e.g., Nakariakov fc Verwichte 2005). 
The maximum difference between the kink speed and 
Alfven speed in going from a uniform density environ- 
ment to a loop in an evacuated environment is a fac- 
tor of a/2- This corresponds to the maximum system- 
atic error in the inferred value of the magnetic field due 
to an interpretation in terms of Alfven waves or MA 
kink mode waves. Recent spectroscopic observations 
with EIS were used to detect a coronal oscillation and 
to map the electron density variation across a coronal 
loop (Van Doorsselaere et alj 120081) . That study indi- 
cated a 20% density enhancement of the loop over its 
surroundings which gives rise to only a 5% difference in 
the value of the inferred magnetic field using either an 
Alfven wave or MA kink mode interpretation. The in- 
vestigation of Young et al.' (2008) into Fe XII and Fe XIII 
electron density diagnostics with EIS indicates discrepan- 
cies of up to 0.5 in the log of the inferred electron density 
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Fig. 8. — Maps of the outward (panel A), inward (panel B) directed and unfiltered (panel C) wave phase speeds and respective errors 
(panels D, E and F). 
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that are due, principally, to uncertainties in the atomic 
mode ls employed by the CHIANTI database (|Dere et al.l 
[1991 . While we acknowledge that it is important to 
determine the actual type of wave responsible for these 
ubiquitous coronal oscillations and to infer physically 
meaningful information from them, we suggest that the 
uncertainty in the magnetic field inferred from observed 
wave phase speeds is dominated by uncertainty in the 
density to a degree that the difference between the kink 
and Alfven speed is small in comparison. In the future, 
CoMP will complement all wave sequences with measure- 
ments of circular polarization (Stokes V) to directly mea- 
su re the LOS magnetic field strengths (see, e.g.. Fig. 8 
of iTomczvk et a0l2QQ8f ). and d ensity diagnos tics of the 
Fe XIII line pair (see, e.g.. lPennlll994l : iSingh. et al.ll2QQ2r ) 

in an effort to a ddress this issue . 

Finally, as in iTomczvk et al.l (|2QQ7l ) , we can estimate 
the energy carried by the ubiquitous coronal waves as: 
Fw = p{'^'^)Vphase^ where p is the density, v is the veloc- 
ity amplitude, and Vphase is the phase speed of the waves. 
Assuming a typical range of electron densities of 10^ — 
10^ cm-3, we find p - 2 • 10"^^ - 2 • IQ-^^g cm"^; also, 
V ~ 300 m/s and the median phase speed over the field is 
549.2 ±28.6 km/s. Then, the fiux of energy propagating 
in the observed waves is Fw = 10-100 erg cm~^s~^ com- 
pared with the estimate of 3 • lO^erg cm~^s~^ required 
to balance the radiative l osses of the quiet solar corona 
(jWithbroe fc N oveslfT977h. We concur with our earlier 
estimate ( Tomcz vk et al.l [20071 ) that the energy carried 
in the resolved ubiquitous coronal waves is insufficient to 
heat the solar corona by at least three to four orders of 
magnitude. Calculation of the energy carried by these 
waves assuming they are MA kink modes could reduce 
the e nerg:y estimate by as much as an order of mag:ni - 
tude (|Van Doorsselaere. Nakariakov fc Verwicht^ l2008f ) 
but leaves our basic conclusion intact. As in our ear- 
lier paper, we caution that this calculation includes only 
the waves that we can spatially resolve with our instru- 
ment and is therefore an underestimate. The ability of 
these waves to heat the corona remains an open question. 

6. CONCLUSION 

Precision Doppler imaging of the corona with the 
CoMP instrument has revealed oscillations which are 
ubiquitous in space and time. In this paper, we have 
exploited the ubiquity of the Alfvenic motions to con- 
struct the first space-time and k — uo diagrams of the 
corona. These allow us to separate outward and in- 
ward propagating waves and obtain a precise estimate of 
the spatial variation of the plane-of-sky projected phase 
speed. The separation of inward and outward waves is 
critical to the precise estimation of the phase speed and 
represents a significant improvement over the analysis of 



ITomczvk et al.l (|2007[ ) . The significant disparity between 
inward and outward wave power over most of the corona 
leads us to conclude that the dissipation time scale for 
these waves is roughly equal to the time of wave propaga- 
tion across the loops. This is consistent with a picture of 
the corona where low-frequency Alfvenic motions suffer 
from significant damping as they propagate along mag- 
netic structures. The power law spectra with a slope of 
-3/2 that we observe is consistent with MHD turbulence. 
However, other dissipati on pr ocesses such as wave mode 
conversion (e.g. Melrosd ll977l ) cannot be ruled out. 

We have shown that observations of this type can pro- 
vide precise phase speed information which is an essential 
ingredient to the broad practical application of coronal 
seismology. However, since we lack information on the 
density of the corona in the region we have studied, we 
can not perform real coronal seismology with this data 
set at this time. We suggest that the lack of precise 
density information in coronal structures is the largest 
obstacle to the widespread application of coronal MHD 
seismology to the study of coronal magnetism. Current 
uncertainties in the density distribution in the corona 
dominate even over uncertainty due to the assumption 
about the nature of the wave being observed. 

Coronal seismology measurements constrain the plane- 
of-sky projected phase speed and therefore the transverse 
component of the magnetic field which is complemen- 
tary to the coronal Zeeman Effect measuremen ts (e.g. 
iLin. Kuhn fc Coulterl [2003 : ITomczvk et al.ll2008[ ) which 
constrain the LOS component of the magnetic field. The 
combination of seismology and Zeeman diagnostics then 
has the potential to constrain the complete vector mag- 
netic field in the corona. We are currently in the process 
of relocating the CoMP instrument from Sunspot, NM to 
Haleakala, HI where we expect the number and quality 
of coronal observations with the CoMP instrument to be 
greatly increased. 
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